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Barriers and Dispersal Surfaces in Minimum-Time Interception^ 


N. RAJAN^ AND M. D. ARDEMA^ 

Abstract. Dispersal points and points on the Barrier are located 
for a class of pursuit-evasion and interception problems. These 
points are determined by constructing cross sections of the iso- 
chrones and hence obtaining the Barrier, Dispersal, and control- 
level surfaces. The method is illustrated through examples from 
minimum-time interception of a target moving in a straight line in a 
horizontal plane. 

Key Words. Barriers, minimum-time interception, aircraft flight 
paths optimization. 

1. Introduction 

The area of optimizing aircraft maneuvers is quite complex and has 
received extensive attention. In Ref. 1, minimum-time turns to a line or a 
point in the plane are investigated when the final velocity and heading are 
specified. The destination point or line is taken far away from the starting 
point so that the aircraft has sufficient time to accelerate to maximum speed 
in all cases. The solution is obtained by integrating the optimal trajectory 
backward from termination. The assumption of long-duration makes the extremals 
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G. Leitmann for their comments. The first author wishes to acknowledge the 
financial support of the National Research Council. 
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3. onG*p3iT3niGt6ir f 3iniXy • TIig two^boundsiry^vsluG piroblsm of choosing the GxtiTG™ 
mal that satisfies the given initial conditions is thus simplified. Under the 
restriction of long—duration maneuvers, time— scale separation of turning and 
acceleration is possible, and singular perturbation approaches (Refs. 2, 3) 
have been successful in obtaining feedback approximations. However, short- 
and intermediate— duration maneuvers do not show up in these analyses. Aircraft 
maneuvers in a tail-chase pursuit-evasion game have been investigated in Ref. 4, 
employing a differential dynamic programming algorithm. The numerical^., 
examples show tail-chase maneuvers of short duration. 

In this paper, minimum— time interception of a target moving in a horizon- 
tal plane is analyzed as a one-player differential game (Ref. 5), The game 
solution gives the interceptor’s strategy in a feedback form. It first 
addresses the question whether the target can be captured from all initial 
positions relative to the interceptor. If not, the Barrier, a closed-surface 
bounding the initial positions from which capture is possible, is constructed. 
Even if capture is possible for all initial states, a Barrier may exist, but in 
this case it is an open surface across which the game value, here the time— to— 
capture, is discontinuous, and the nature of the optimal trajectories changes 
drastically. 

The game solution maps the controls as a function of the state within the 
capture region. This is done by backward integration of the state and adjoint 
equations from termination. The integration has to be stopped at surfaces 
where the ad joints are discontinuous; such surfaces are called singular surfaces 
(Ref. 5). One of the most frequently encountered singular surfaces is the 
Dispersal surface. Along such a surface, one player can choose between two 
distinct strategies and still obtain the same payoff. Constructing the .Barrier 

and Dispersal surfaces is therefore an inherent part of a differential game 
solution. 
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In a game with four or five state variables, visualization of the' Barrier 
and Dispersal surfaces is difficult. This is circumvented by constructing the 
cross sections of these surfaces with two or more of the state variables held 
constant. A stack of such cross sections then describes the surfaces. This 
paper describes a method of computing the locations of points on these cross 
sections. Although minimum-time interception in the plane is used as an 
example, the extension of the method to pursuit-evasion in the plane, which 
is straightforward, is also outlined. 

The analysis uses a device first developed in Refs. 6-8. The motion of 
the interceptor and the target are described in a coordinate system whose 
origin is fixed at the interceptor’s terminal position, and the x-axis is 
along the terminal line of sight. The equations that describe' the inter- 
ceptor's optimal motion backward in time from termination do not contain any 
reference to the target motion. They can be integrated independently to give 
extremal trajectory maps (ETM) (Ref. 6). Using this concept, cross sections of 
the Barrier were constructed by direct iteration of the terminal Mach numbers 
of the pursuer and evader in Refs. 6-8. In this paper, constant initial Mach 
number and relative heading cross sections of the isochrones (constant-time 
loci) are drawn. At Dispersal points (DP) , two branches of an isochrone 
intersect. An isochrone ends if capture becomes impossible for the given time- 
to-go. The endpoints thus fall on a Barrier cross section. 

The interception problem is formulated and the necessary conditions for 
optimality are derived in Section 2. Extremal trajectory maps are shown for a 
given initial interceptor speed. The method of constructing a global solution 
is described in Section 3. Isochrone cross sections for interception are 
depicted, and some interception encounters in the plane are shown. The exten- 
sion of the method to pursuit-evasion is outlined. 
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2 . 


Problem Statement 


2.1 Formulation. The motion of the target in the plane is described in 
terms of its coordinates [x^(t), y^(t)] in any conventient sec of axes. 

Since the motion is prescribed a priori, the target motion is a function of 
time alone. In the same coordinate frame, the Interceptor's motion is 
goverened by 

Xj = MjCosBj, (1) 

y^ = MjSinBj, (2) 

Mj = A(Mj)tt - B(M^) - (3) 

and 

Bj = o)a(M^)/Mj, (4) 

subject to the state-variable constraint, 

M^e[M, H], ( 5 ) 

with the throttle setting tt and the bank control o) chosen such that 

t£[ 0, 1] and w£[-l,l]. 

The above description of the aircraft's motion is the same as that given 
in Ref. 6. Briefly, Bj is the interceptor's heading measured from the x-axis 
(counterclockwise positive); is the Mach number; A represents the maximum 
thrust, B is the zero bank drag, C is the lift-induced drag at the maximum 
bank, and a is the angular acceleration. The target is considered to have 
been intercepted when brought within a line of sight (LOS) distance of R of 
the interceptor. This is. 
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( 6 ) 


[x^(tp -R^ 

at the instant of interception t^, the interceptor must be reducing the 
separation distance or 

r = M^(tj)cos6^(t^) - Mj(t^)cos6j(tj)<0 (7) 

to effect capture. Where interception occurs, the payoff is the time taken to 
capture the target. The performance index, augmented by the terminal con- 
straint (6), is 

t 

= v[(x^(t^) - Xj(t^))2 + (y^(t^) - y^(t^))^ - R^] +/ ^dt. (8) 

As in the analyses of pursuit-evasion, (Refs. 6-8) the coordinate axes, 
to which the above equations refer, are selected such that the origin coincides 
with the interceptor's terminal position, and the x-axis is along the terminal 
line of sight (Fig. 1) . The position and orientation of the axes thus change 
from engagement to engagement. As is shown below, this choice of axes makes 
the trajectory equations and the expressions for the optimal controls for the 
interceptor independent of the target motion. The analysis is simplified and 
leads very naturally to the method of locating Dispersal points, which is 
developed later. 

2.2 The Necessary Conditions. The Hamiltonian for the minimum- time 
interception problem is written as 

H = 1 + X^x^ + X^y^ + (X^ + u)M^ + XgBj., (9) 

where X^, X^, are the position adjoints; X^ and Xg are the velocity and 
heading adjoints; and the multiplier y accounts for the first-order state- 
variable (Ref. 9) constraint on the interceptor's Mach number. It satisfies 
the Kuhn- Tucker conditions. 
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<0 if =M, 

li = 0 if M<M^<ir, (10) 

>0 if = H. 

From the Euler-Lagrange equations and the transversality conditions plus 
the particular choice of coordinate axes mentioned above, it can be shown that 
the ad joints satisfy (Appendix) 


X 

X 


-2vR, 


( 11 ) 


X 

y 


EO, 


( 12 ) 


and 


Xg = -2vRy^. ( 13 ) 

The vanishing of the Hamiltonian at t^ gives the value of the multiplier v 


V = (2R(M^^cos6j^ - M^^cosS^j))>4; (14) 

V is positive because the target's velocity relative to the interceptor has a 
negative component along the line of sight at the endpoint. 

Because H does not involve time explicitly, a first integral of motion 
is given by (Appendix) 


min H = -2vR M ,cos6„-. (15) 

(•^»w) 

When the normalized velocity adjoint and the multiplier y are defined, 

repectively, as 


A 


(16) 
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and 


ii = u/2vR, (17) 

and by using (11)-(13) in (9), then (15) becomes 

min H •= -M ,cos6_,, (18) 

(Tr,o))^ “ “ 

where 

Hi ^ (p^+u)Mi-Xi-yi®I. 

Minimizing with respect to subject to the constraints '^"£[0,1], 

0 ) 6 :[-l,l] yields the following expression for the optimal controls: 

7T = 1 

0 ) = sat(-ay^/(2CM^(p^+y))) for P^+y<0 (20) 

and 


7T = 0 

03 = sign(y^) for Pj^+U>0. (21) 

The normalized velocity adjoint appearing in the above expressions satisfies 
the differential equation 

pj^ = cos6i-(p^+jj) (-dA/dMi-oj^dc/dMi)+ (22) 

0)yi(da/dMi-a/Mi)/Mi. 

Note that the canonical equations (1), -(A), and (22), together with 
the controls (20) and (21) contain no reference to the target motion; the same 
type of "open-loop optimality" found in pursuit-evasion is obtained here. 
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Backward integration of the extremals depends only on the terminal-velocity 
vector of the interceptor and is independent of target motion. Thus a map of 
the extremals may be generated integrating backward from the origin; this is 
called the extremal trajectory map (ETM) in Refs. 6 and 7. 

2.3 Extremal Trajectory Maps. Each extremal in an ETM may be consid- 
ered to be specified by a triplet of values (6^, M^, t^), where 6^ and M^ 
are the terminal heading and Mach number, and x^(= t^-g) is the time-to-go. 

For any given pair of values for and it is possible to select an 

M^ value, through a one-dimensional search, such that the Mach number at x^ 
matches a given value M^. The ETMs used here are collections of extremals 
that start out from a given initial Mach number. They can be combined with a 
target/pursuer/evader's ETM to yield a global solution for any given 
encounter. For this, apart from the triplet (g^, M^, x^) , the initial 
coordinates (x,y) and heading g are needed. The behavior of the bank and 
throttle along the extremals determines the values of x^ for which iso- 
chrones need be drawn (Section 3) . 

The trajectories in the plane for = 10 sec are shown in Fig. 2a. 

The variation of Mach number, heading and bank with retro time is shown in 
Figs. 2b - 2d. For g^ = 0°, the aircraft flies in a straight line at full 
throttle, accelerating from Mach 0.9 to Mach 1.02. The extremal with 
Sj = 0.33® starts out at half-bank and gradually straightens out to near- 
level flight. The angle of turn Ag = 23®, and the terminal Mach number is 
slightly above unity. On the next extremal, the bank decreases from unit 
magnitude at start to below 0.1 at termination; Ag = 45®, and the average 
turn rate is about half the maximum instantaneous at Mach 0.9. As g^ 
increases beyond 0.85®, the bank saturation instant moves closer to termination, 
till at g^ = 15®, the entire trajectory is flown at full bank. 
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These trajectories represent different trade-offs between turn and 
speed. However, their initial points would not appear on the cross section of 
the unit bank control surface taken with the initial Mach number equal to 
0.9. Between “ 15® and 77®, the extremals are constant control arcs; both 
bank and throttle are full throughout. The Mach number time-history is 
identical for all extremals in this range of 6^ and is hence a single curve 
in Fig. 2b. The turn angle A6 » 109®, or the average turn rate = 10.9®/sec. 
Between 77® and 90®, the extremals start out with the throttle closed and 
switch to full throttle sometime later; this switching instant moves toward 
with increasing Sj till, at 90®, the extremals terminate with zero 
throttle. The region between ~ 90® and 100®, shows throttle switches from 
zero to full and back to zero; .for 6^ = 100°, at T = 9.5(0 to 1) and 
6.6(1 to 0) sec, respectively. Beyond = 100°, the extremals are all 
flown at zero throttle, full bank; the Mach number variation is independent of 
6^ (Fig. 2b). 

The 140° and 160® extremals show switchings in bank (Fig. 2c) where 
they intersect the x-axis. Such maneuvers reduce the turn angle; the net 
effect is to brake down while heading in the same direction. Such a maneuver 
is resorted to when the target is close behind and heading in the same direc- 
tion. In Fig. 2a, the extremals for 6^ * 140® and 160® appear to intersect. 
However, the heading angles at the point of intersection are different, as can 
be seen from Fig. 2c. 

Bank and throttle switching disappear as increases beyond 25 sec. 

The behavior of the extremals with increasing becomes more uniform and can 

be summarized as follows: as x^ increases, the values of B^ for the same 
AB decrease steadily (for x^ = 60 sec, B^ = 0.0005® for AB = 200°); all 
extremals have partial bank segments and are full— throttle throughout for the 
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initial Mach number of 0.9. As increases further, the terminal Mach 

number hits the speed boundary. The extremals then consist of an uncon- 
strained (maximum Mach number constraint) arc accelerating up to M, a 
circular arc at M, then a straight dash at M along the negative x-axis. 

Such extremals, called cruise arcs, are similar to those discussed in Ref. 1. 
However, the use of these extremals to locate Dispersal points depends only on 
a decoupling coordinate transformation approach being possible; it is discussed 
next . I 

3. Construction of a Feedback Solution 

A feedback solution to interception or pursuit-evasion must give the 
control strategies and the value at any point in the state-space. One way of 
the solution is to draw the boundaries of the regions within which 
the strategies are constant. Where a control varies continuously over a range, 
as does the bank here, the surfaces at which the control attains specified 
values are mapped. Since these boundaries and surfaces are in a four or five 
dimensional state-space, they can be visualized by cross sectioning. They are 
sectioned by keeping the initial Mach numbers of the interceptor and target, 
or of the pursuer and evader, and the initial relative heading, constant. The 
cross sections are then plotted relative to the interceptor/pursuer and become 
curves in the plane of the encounter. 

The control-level surfaces. Barrier Points (BP), and Dispersal Points 
are all located by drawing cross sections of the isochrones (the optimal parts 
of the isochrones are isocost). The control-level surface links the points on 
the isochrone at which the controls take on a specified set of values. 

Dispersal Points are located at the intersections of different branches of the 
same isochrone. Where the isochrone cross sections are not closed, their 
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endpoints are on a Barrier cross section. Such a section separates two 
families of isochrones. The construction of the isochrone sections for inter- 
ception is described next. 

3.1 Isochrones for Interception. Assume that the interceptor and 
target ETMs- are as shown in Figs. 3a and 3b, respectively. The interceptor's 
map comprises extremals generated for a specified time to go with the terminal 
heading 6^ as a parameter. For each extremal in the map, the initial Mach 
number is the initial heading measured relative to the terminal line of 

sight equals The target's initial Mach number is given as Since 

the target's motion is prescribed, the target's terminal Mach number M^^ 
and the angle of turn at are known. Although the target's motion 

in the plane is fixed, its heading relative to the terminal line of sight 3^^ 
varies with the different paths flown by the interceptor. The target's ETM 
consists of the target's path in the plane rotated through several different 
values of (the term ETM is used only as a convenience). 

If the relative heading of the isochrone section being constructed is 
then the initial heading of the target measured from the terminal LOS must 
be equal to Since the target's angle of turn is AB^, its terminal 

heading equals 3j^+3^-A3^. In the construction of a locus point (Fig. 3c) , 
the interceptor's extremal is first laid off as lA. The terminal LOS is AB. 
Since the target's heading relative to AB is known, it can be laid off as 
shown, TB. The point T is the point on the locus, provided the terminal LOS 
rate is negative; that is 

M^^cos(B3.Q+6Q-A3^)-M^^cos3f<0. 

Different values of 3^ give other points on the locus. Points of 
interest on the extremals, for example bank saturation, throttle, and bank 
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switching, also get mapped on the isochrone locus section. For a section with 
a different relative heading, these points are rotated by the difference in 
method plots all candidate points for the isochrone. This is true 
because extremals for all possible values of 6^ are generated. Some of the 
resulting points are not optimal because there are other constant~time loci 
points for a lower time-to-go overlapping them. A Dispersal point occurs if 
two extremals with different values of 3^ pass through the same point in 
reduced space. Neither extremal is globally optimal beyond the DP. 

The location of Dispersal and Barrier points and the construction of 
isochrones are explained by examples. 

3.2 Examples. Figure 4a shows some isochrone section for an inter- 
ceptor initial Mach number of 0.9 and a relative heading of 135°. The points 
^1 ^4 those for which 6^ = 0° ; the interceptor just flies straight 

ahead at full throttle to meet the target. To the right of these points, the 
interceptor banks to the right, and the terminal heading 6^ is positive. 

For all the sections, 6^ increases in the direction C to B. The sections 
do not extend beyond the points B and C because the "reachability” 
condition is not met; thus, the envelope of these points delineates the 
3 q = 135°, M^^ = 0.9 section of the Barrier (Ref. 5). However, since the 
interceptor has a maximum Mach number capability of 1.6, it can capture the 
target, given sufficient time. This can be seen from Fig 4b where the 30 sec 
locus is also shown. Across B^, there is a discontinuity in the time-to-go. 
From B^ to P, the time— to— go is 30 sec. The PQD^ portion of the constant- 
time locus for t^ == 30 sec is nonoptimal and does not form part of any iso- 
chrone. All points on PQD^ are conjugate points; however, only D^ and P 
actually appear on the isochrone sections. This suggests that not all values of 
0^ lead to extremals that are globally optimal for a given value of 6^; 
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there are points for = 20 sec which fall between the origin and the 15- 
sec locus and are hence obviously non-op timal globally. 

Figure 4c shows a more complete picture of the isochrone sections with 
ranging from 5 to 60 sec. The 60-sec locus is the first one that is 
closed. Even for an initial range of 2 km, 60 sec is the globally optimal 
time-to-capture because of the unfavorable initial conditions. The target 
has a higher initial speed and its direction is away from the interceptor. 

The interceptor turns left and accelerates to Mach 1.2 with its terminal 
heading almost aligned with the final LOS. This encounter starting from P 
is depicted in real space in Fig. 5. The locus shows a ^Dispersal point at 
D. Extremals with two different values of 6^ pass through D; they are 
shown in the plane in Fig. 6. 

In Fig. 7, an isochrone section with = 0® and = 20 sec is 
shown. The Dispersal point at is one along a Dispersal surface where the 

interceptor flies symmetric paths starting with either a left or a right turn; 
the values of at D^^ are +116®. The encounter at D^^ is shown in Fig. 8. 

Another interesting kind of Dispersal point occurs at D 2 in Fig. 7. .Two 
branches of the locus, one with 3^ decreasing toward 90® and the other with 
3^ increasing, intersect at D^. The encounter at D 2 is shown in Fig. 9. 

The additional computational effort to map the feedback solution, given 
the ETMs, is relatively small. Moreover the logic of determining the Barrier 
and Dispersal Points can be automated. The generation of the ETM data storage 
is a straightforward, though tedious process; once done, however, the data are 
common to any interception/pursuit evasion problem in which the given aircraft 
is a participant. The method of constructing maps of strategies and isochrones 
can serve as a useful tool for design, displays, training, and testing of sub- 
optimal approximations. The extension of the method to interception and 
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pursuit-evasion in three dimensions is, in principle, straightforward. If 
the aircraft is modeled based on the energy-state approximation, with the 
specific energy replacing the Mach number as a variable, a similar structure 
is obtained. The extension of the method to pursuit-evasion is outlined. 

3.3 Pursuit-Evasion. The above method of constructing isochrone 
sections can be employed for an actively evading target. The cross sectioning 
would be done for fixed initial Mach numbers and relative heading for the 
pursuer and evader. The evader's ETM would be constructed in exactly the same 
way as the pursuer’s, with 6^, T^, as parameters and M^ iterated till the 
Initial Mach number is matched. For each point on a constant-time locus, the 
evader’s ETM would be searched for the extremal that has an initial heading 
equal to The "reachability” condition would be checked for the 

extremal. This search for the extremal in the evader's ETM is independent of 
the search for the M^ value that generates it. Other features of the method 
are the same as for interception. 

4. Conclusions 

A method of locating Barrier and Dispersal points for a class of minimum- 
time interception and pursuit-evasion problems has been developed. The method 
obtains these points, which are an essential part of the game solution, by 
drawing cross sections of the isochrones. Examples are shown for minimum-time 
interception in a horizontal plane. The method can be extended to solve pursuit 
evasion between two aircraft in the plane and in three dimensions. 
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Appendix 


The differential equations for the ad joints are 

\ = 0, (Al) 

- 0, (A2) 

* -X^cos6j-XySin8j-(Xj^+y)(TdA/dM-dB/dM-a) ^dC/dM) (A3) 

+i.a3(a/M-da/dM) /M, 

P 

and 

^3 = VrVi. 

For the choice of coordinate axes shown in Fig. 1, the transversality 
conditions give 


X^(tj) = -2vR, (A5) 

Xy(tf) = 

and 

H(tj) = -2vRHj,fCOs6^f. (A7) 

In (A7) , and 3,p£ are the target's terminal Mach number and heading. 
Equations (Al) , (A2) , and (A4) , together with the end conditions (A5) and (A6) , 
yield (11) - (13). The substitution of (A5) in (A7) leads to (14). 
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Fig. 5 
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SHOWN IN REAL SPACE 
POINTS MARKED AT 10 sec INTERVALS 

□ INTERCEPTOR PATH- /3^ = -4.1 e- 4® M^ = 1.209 
O TARGET PATH 135° Tf=60sec 
A INTERCEPTOR PATH - = 4.3 0 - 4° Mf = 1.21 


10 km 



Fig. 6 




SHOWN IN REAL SPACE 
POINTS MARKED AT 1 sec INTERVALS 

□ INTERCEPTOR PATH - = 1 16° Mf = 0.576 

O TARGET PATH - /3q = 0° = 20 sec 



Fig. 8 
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SHOWN IN REAL SPACE 
POINTS MARKED AT 1 sec INTERVALS 

□ INTERCEPTOR PATH - = 99.5® Mf = 0.576 

O TARGET PATH - = 0® tf = 20 sec 

A INTERCEPTOR PATH - /3f = 83.5® Mf*0.8 


1 km 
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